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^"v We have made an analysis of limits on the neutrino mass based upon the formation of large-scale 

structure in the presence of a primordial magnetic field. We find that a new upper bound on the 
neutrino mass is possible based upon fits to the cosmic microwave background and matter power 
spectrum when the existing independent constraints on the matter density fluctuation parameter 
us and the primordial magnetic field are taken into account. 
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o 

J3 '■ I. INTRODUCTION 

The evolution of the cosmic matter density field and large-scale structure (LSS) can be affected by a finite neutrino 
+3 ', mass [l|. For example, when the velocities of finite-mass neutrinos are large, the growth of density fluctuations will 
^3 ' be impeded [2| on the scale of neutrino free-streaming. Cosmological constraints on the neutrino mass are thus of 
considerable interest. Recently an upper limit to the neutrino mass of order < 0.1 — 1 eV has been deduced from a 
combination of tritium beta-decay endpoint experiments | 3| a nd cosmological observations [4|-|7| • 

At the same time, magnetic fields have been observed [cHLlj in clusters of galaxies with a strength of 0.1 — 1.0 [i, G. 

One possible explanation for such magnetic fields in galactic clusters is the existence of a primordial magnetic field 

(PMF) of order 1 nG whose field lines collapsed as the cluster formed. A varietyof mechanisms for the origin and 

\^q \ amplification of a PMF on the scale of galaxy clusters havs been proposed |12h14| and the associated PMF could 

in 



have influenced a variety of phenomena in the early Universe 15] such as the cosmic microwave background (CMB) 
|16l - l39| , and the formation of LSS [40-45J . Of particular relevance to the present work is the fact that the effects of a 
finite neutrino mass on the matter density fluctuations are degenerate with the effects of a PMF. Therefore, these two 
effects must be constrained together when using the fluctuations of the matter density to find limits on the neutrino 
mass and the PMF parameters. 

In this regard, the alternative normalization parameter as is of particular interest as a measure of large-scale 

structure effects. It is defined [46| as the root mean square of the matter density fluctuations in a comoving sphere 

JC ■ of radius 8/i _1 Mpc. It is determined by a weighted integral of the matter power spectrum. Observations which 

determine as provide information about the physical processes affecting the evolution of density-field fluctuations and 

the formation of structure on cosmological scales. 

The detailed mechanisms by which a PMF can affect the density field fluctuations on cosmological scales has been 



described in Refs.[42t |44|. Of course, as is also affected by the presence of a PMF. Indeed, in [45j we demonstrated 
for the first time that there was a degeneracy between the parameters of the PMF and the neutrino mass. 

In this article, we expand on our previous work 45] by also considering the influence on the matter contributions 
|38l |39| . By considering the effects of a PMF and finite- mass neutrinos on as and comparing theoretically deduced 
values for as with the observed range, we obtain not only insight into the underlying physical processes of density 
field fluctuations in the presence of a PMF, but also show that the constraint on the sum of neutrino masses from as 
changes when the effects of a PMF are included. 
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II. THE MODEL 

To evolve the primary density perturbations in the presence of a PMF we begin with magnetized isocurvature initial 
conditions for the fluids. For the present purposes we also fix the cosmological parameters to those of the best-fit flat 



ACDM model as given in the WMAP 5 yr analysis in Ref. |47|, i.e., fl b = 0.0462, tl c = 0.233, n s = 0.96, h = 0.701 
and r c = 0.084, where fif, and 51 c are the baryon and cold dark matter densities in units of the critical density, ns is 
the spectral index of the primordial scalar fluctuations, h is the Hubble parameter in units of km s Mpc , and r c 
is the optical depth for Compton scattering. 

The ionized baryons and electrons arc influenced by Lorentz forces when a PMF is present. Photons are then 
indirectly affected by the PMF through Thomson scattering before the epoch of photon last scattering. We assume 
that the PMF was generated earlier during the radiation-dominated epoch. We use a Friedmann-Robertson- Walker 
(FRW) background cosmology for the linear perturbations. For the time evolution we treat the energy density of the 
PMF as a first-order perturbation and a stiff source. Therefore, all back reactions from the fluid can be discarded. 
We also assume that the conductivity of the primordial plasma is very large and that the electric field is negligible, 
i.e. E ~ 0. This "frozen-in" condition is a very good approximation (18j ] . On very large scales the time evolution of the 
PMF can be decoupled from its spatial dependence, i.e., B(ry, x) = B (x)/a 2 . This leads to the following simplified 
electromagnetic energy-momentum tensor: 



.oo ,„ ., _ W 



T°Wx,t) = ^- , (1) 

T lO [E M](x,T)=T ofe [E M](x,T)=0 , (2) 

-T* [EM] (x,r) = erf = J- {2B\ X )B k (x) - <5 lfe i?(x) 2 } , (3) 

where we use natural units c = h = 1. For a PMF that is statistically homo gene ous, isotropic and random, the 
fluctuation power spectrum can be parametrized as a power-law P(k) oc fc" B [18L |28| where «b is the spectral index. 
A two-point correlation function for the PMF can be defined [l8| by 

. (J \n B +8 f>2 

B*(k)BJ (fc')) = \ n+3 ^ r * k n *pv(k)5(k-k'), k<k c , (4) 



2k n + 3 p ( n B+3 ^< 

where B\ = |Ba| is the strength of the magnetic comoving mean-field derived by smoothing over a Gaussian sphere 
of comoving radius A and k\ = 2ir/A (with A = 1 Mpc in this paper), and P tJ (fc) = S 1 ^ — k % W /k 2 . There is a natural 
cutoff wave number kc in the magnetic power spectrum which is dependent upon the damping scale of the PMF by 
radiative viscosity. It is defined in Refs. [48j-l50J. Hereafter, we work in k-space and denote all quantities by their 
Fourier transform convention F(k) = J d 3 xex-p(ik ■ x)F(x). 

We evaluate the PMF source power spectrum using the numerical methods described in Refs. [32|, |42[. Using this, 
we can quantitatively evaluate the time evolution of the cut off scale and thereby more reliably calculate the effects 
of a PMF. 

We adopt the model of Ref. [33, [38| for estimating the effects of finite- mass neutrinos and the PMF on fluctuations of 
the matter density fields in the early Universe. We also utilize adiabatic initial conditions for the matter contributions 
as in Ref. [38, [39| . This leads to stable numerical calculations of the curvature perturbation of the scalar mode for 
all scales and times. This is an improvement over previous numerical estimations for which the scalar curvature 
perturbations were too small to stabilize numerical calculations for large scales and early times. Thus, we can obtain 
consistent results for all scales and times of interest in the present work. 

III. RESULTS AND DISCUSSIONS 

Details of the effects of a PMF on the cosmological density field fluctuations can be found in Ref. [42j . Here we 
briefly review these effects and show that a relation exists between erg and the PMF parameters. The observational 
constraints [5ll - [54j | on erg, therefore, lead to a strong constraint on the PMF parameters and the neutrino mass. Since 
the parameters of the PMF have a strong degeneracy, the existence of a Bayesian prior constraint on eg can be used 
to effectively constrain the PMF. Also, since eg is constrained by diverse observational data on linear cosmological 
scales, we can obtain a reliable prior for use in determining the probability distributions for the parameters of the 
PMF from CMB observations. 

The PMF effects dominate the matter power spectrum for wave numbers k > O.l/i Mpc -1 [42] (Fig. 1). This is because 
the PMF energy density fluctuations depend only on the scale factor a and can survive below the Silk damping scale. 



Therefore, the PMF continues to be a source for the fluctuations via the Lorentz force even below the Silk damping 
scale. In the case of no correlation between the PMF and the matter density fluctuations, the matter power spectrum 
is increased by the PMF independently of whether the PMF is dominated by the pressure or the magnetic tension [42j. 
A PMF affects the power spectrum function P(k) and the matter density fluctuation 5 differently. The total density 
fluctuation 6 can be smaller or larger depending on whether the effect of the PMF is mainly due to its pressure or 
tension. On the other hand, the power spectrum function P{k) always increases in the case that the PMF does not 
correlate with the primordial density fluctuations. This is because P(k) oc 8 2 and is not affected by the sign of 5. 

A. PMF parameters and ag 

The alternative normalization parameter <r% can be determined from a study of the physical processes of density 
field fluctuations on cosmological scales within the linear regime. Recently as has been constrained in this way by 
observations [5lT - l54| to be in the range 0.7 < as < 0.9. From this we can obtain strong constraints on the PMF 
parameters by numerically calculating a$ under the influence of a PMF. 

We have shown [29|, |32|, |42[ that the possible discrepancy between theoretical estimates and observational temper- 
ature fluctuations of the CMB for higher multipoles (I > 1000) could be attributed to a PMF of comoving strength 
1.0 nG< B x <2.0nG. Since other constraints on the amplitude of the PMF imply B x <10 nG flU-Hl I2U [30. [55Jj . our 
allowed range is consistent with previous works. 

The value of a% derived by a CMB analysis including such a field strength for the PMF is also 0.78 < a% < 0.84. 
As noted above, in this work we consider a% constraints from both fits to the CMB and from large-scale structure. In 
this case, we adopt a prior for as in the range of 0.75< as <0.85 as a reasonable average of constraints from observed 
structure and fits to the CMB. 

B. Constraints on the neutrino mass and PMF parameters from as 

Figure [2] shows the constraints on the PMF parameter B\ and the sum of neutrino masses Em„ for various fixed 
values of ub and ranges of as as labeled. Since the PMF power spectrum depends on ub, PMF effects on the density 
fluctuations for small scales decrease with lower values for tib- Based upon previous work ka-u( the upper limit on 
the total neutrino mass is expected to be of order < 0.1 — 1 eV. For such masses the neutrino velocities can be quite 
large while the mass density of neutrinos is still a significant fraction of the total dark matter density. Therefore, the 
growth of density fluctuations on the free-streaming scale of neutrinos will be hindered. 

Thus, neutrinos decrease the matter density fluctuations, while a PMF increases the matter density fluctuations 
(cf. Figure [T]). Furthermore, for Ub within the ranges constrained by previous work [20, (2y|, a PMF with B\ >1 
nG affects the matter density fluctuations (see Fig. [5J . When the mass of the neutrinos is constrained from matter 
density fluctuations in the presence of a PMF, the deduced upper limit is heavier with a PMF than without. 

The expected parameters of the PMF from the combined analysis of the CMB and observed magnetic fields in 
galactic clusters is B\ <2.0nG(la) and <3.0nG(2cr) [23, 123, |37[ , while the expected value of as based upon observations 
is 0.75< as <0.85 as noted above. For this range of as, the sum of the neutrino masses is constrained to be 
J2n m v < 0.8 eV (from the la constraint on B\) and < 2.2eV (from the 2a constraint on B\) for N v = 3. This is a 
larger upper limit than that deduced previously because the effect of the PMF cancels the effect of neutrinos on the 
density fluctuations 0-0] ■ 

C. Future Work 

Komatsu et al. (2008) noted that the effect of a finite neutrino mass is very small for the detection range of the 
WMAP, e.g. £ < 1000 [6]. Therefore, as a first step toward constraining the neutrino mass together with the PMF 
we can use the WMAP CMB power spectrum for t < 1000 to obtain values for the primary cosmological parameters 
without concern for a degeneracy between these parameters and the neutrino mass. One cannot constrain the mass of 
neutrinos very well from such CMB data alone. Hence, one must use other observational data which arc independent 
of the CMB analysis to constrain the neutrino mass. For this purpose, as is very useful since it is determined from 
independent observations of LSS in the regime of linear perturbation theory. Since the main purpose of the present 
work is to examine how the upper limit of the neutrino mass determined from as changes when a PMF is included, 
we have focused on a theoretical analyses of effects of both the PMF and the neutrino mass on erg and the matter 
power spectrum. However, theoretical estimates of as depend upon the cosmological model employed and as has a 
significant degeneracy with Q m = flf, + flcuu, ns, and As, even if such primary parameters are well constrained 



by the CMB data, e.g. WMAPtQ]. Furthermore, there is the degeneracy between the mass of neutrinos and Q m 
[56| . Therefore, one should ultimately consider the degeneracy between the PMF and other cosmological parameters 
including the mass of neutrinos. In subsequent work we will carry out a full-constraint on the mass of neutrinos using 
a Markov chain Monte Carlo (MCMC) method applied to both the CMB and the LSS data. In that effort we will 
need to consider the model dependence and degeneracy between cosmological parameters including the neutrino mass 
and PMF parameters. 

IV. SUMMARY 

Previous numerical estimates of curvature perturbations without matter contributions were too small at early times 
to stabilize on large scales. In the present work we have stabilized the time evolution of curvature perturbations 
in the presence of a PMF while considering the matter contributions. We find that a PMF significantly affects the 
constraint on the sum of neutrino masses. We confirm that the upper limit on the neutrino mass from as in the 
presence of a PMF is heavier than without a PMF even if we consider the matter contributions. We also have shown 
that the prior limited range on the sum of neutrino masses and PMF parameters is within the expected range for eg 
from observations of the LSS. Since the sum of the neutrino masses has a degeneracy with some of the cosmological 
parameters, e.g. f2 m , one should ultimately utilize a more complete MCMC method with the CMB and LSS data to 
resolve such degeneracies and obtain a more accurate constraint on the sum of the neutrino masses. Constraints on 
the PMF parameters and as are soon to be improved by future cosmological observational programs such as Quiet, 
Planck, and the Sloan Digital Sky Survey (SDSS). In principle, by applying our method to such data it will be possible 
to obtain not only the upper but also the lower limits to the neutrino mass from cosmology in the presence of a PMF. 
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FIG. 1: Comparison of the effects of a PMF and finite-mass neutrinos on the matter power spectrum. The left panel is for 
the full power spectrum on scales O.OlMpc -1 < k/h < 0.2Mpc _1 and various parameters as labeled. The right panel shows 
the expanded spectrum in the range of O.lMpc -1 < k/h < 0.2Mpc -1 . Curves and dots with error bars in all panels are as 
indicated in the legend box. In order to better illustrate the effects of a PMF and a finite neutrino mass on the matter power 
spectrum, we have utilized a larger value for the PMF amplitude B\. This figure shows that the matter power spectrum for 
k/h > O.lMpc -1 is affected strongly by a PMF. The magnitude of the effect depends on both the PMF amplitude B\ and 
the power spectral index hb . On the other hand, the mass of the neutrinos decreases the total amplitude of the matter power 
spectrum. Also, the effects of a neutrino mass for smaller scales (large wave number k) are greater than for larger scales (small 
fc). 
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FIG. 2: Excluded and allowed regions in the parameter plane of PMF amplitude B\ vs mass of neutrinos ^2 N _ 3 m v . Shaded 
regions in all panels are allowed by the indicated ranges for as and B\ . Various panels are for the indicated values of the PMF 
power-law spectral index ub- 
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